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On the basis of the perturbation expansion for electron gas energy in the third order 
to the electron-ion potential the pair and irreducible three-ion interaction potentials in 
metallic hydrogen are calculated. The irreducible potential of three-ion interaction has at- 
tractive nature at short interionic separation and oscillates at large ones. The anisotropic 
character of the three-ion interaction is shown. The potential relief of the ions pair relative 
to the third ion is constructed. This relief has some potential wells and valleys which con- 
nect them. The important role of the irreducible three-ion interaction in formation of the 
local order in three- and four-ion clusters is shown. The potential relief of an equiangular 
ion triplet relative to the fourth ion is calculated. This relief has deep potential well for 
fourth ion at interionic separation corresponding interproton separation in molecule H2. 
The quasiclassical probability of ion transition into this well is evaluated. The life time 
of metallic phase of hydrogen relative to the tunneling nucleation of the H2 molecules is 
estimated. 
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1. Introduction 

The significance of manyparticle interactions and their role in forming the struc- 
ture and properties of various condensed systems are not well known yet [|IJ. For 
a calculation of manyparticle interaction potentials the knowledge of the linear and 
nonlinear response functions of electron gas is necessary. At present the explicit 
form of three-pole diagram expression in ring approximation for homogeneous elec- 
tron gas is known 0, @. It was possible to calculate the equilibrium structures 
and dynamic properties of simple metal crystals and metallic hydrogen (MH) in the 
third order to the electron-ion potential (EIP) 0-0, to study three-ion interaction 
potentials in simple metals 

The study of the possibility of metastable state of MH at zero pressure is of 
special importance |4j]. In a case of solid MH calculations have been performed 
by Brovman and Kagan. They calculated the structure, elastic properties, phonon 
spectrum of the solid MH at zero pressure and proved the local stability of the 
metastable MH phase in framework of the manyparticle theory of metals @]. As it 
turned out in the third order to the EIP the energy minimum exists for twoparameter 
hexagonal lattices with triangular string structure Q . For all structures with energy 
minimum the elementary cell volume was Qq = 20, 8a% (clb is the Bohr radius ) and 
interproton spacing fixed along Z axis was equal to d — 2, 04a£. The investigations 
of the higher order approximations did not change conclusions essentially. It is 
connected with the fact that eventually the small parameter in the perturbation 
expansion for MH energy at P=0 approximately equals to 1/5 H,|§. In such a case 
the sum of all perturbation expansion terms with orders n > 5 gives an error as the 
dielectric permittivity of the uniform electron gas does || , || . 
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The estimate of the life time of the metastable MH phase relative to spontaneous 
quantum tunneling transition in the insulating phase is among the most important 
problems. This estimate cannot be made without consideration of the specific mi- 
croscopic mechanism of the nucleation of H2 molecules or ions in metallic phase. 
For the MH as a system in which the tendency to diatomically ordering exists ||, 
the calculations of interionic interactions potentials and elucidation of this tendency 
on the language of groups ion interaction potentials in configurational space have 
great interest. In the quantum kinetics of the new phase formation two characteris- 
tic times may be distinguished In the first place the time of transition through 



potential barrier (in our case the life time of the metallic phase relative to the nu- 
cleation of the H2 molecule). In the second place the total time of the new phase 
formation. As a consequence of a significant difference between electron and proton 
masses during the process of molecular hydrogen phase homogeneous nucleation the 
characteristic time will be determined by slow ionic motion in potential relief created 
by nearest neighbouring ions group. Diatomically ordering in MH corresponds to 
the possibility of the drawing together of ions to the distance corresponding to the 
interatomic separation in molecule H 2 . Inasmuch the direct interionic interaction 
corresponds to the Coulomb repulsion the possibility of drawing ions together may 
be ensure only by means of indirect ions interaction through surrounding electron 
gas. 

The purpose of this paper is to consider pair- and three-ion interactions in 
metastable MH at P=0. The consideration is based on the manyparticle Brov- 
man and Kagan approach 0-0. The peculiarities of the potential relief for ions 
group are discovered and interpreted from the point of view manyparticle tunneling 
formation of the molecular hydrogen phase. 

In sec. 2 expressions for many-particle potentials are received. In sec. 3 and sec. 
4 results of calculations are shown. 

2. Many-body interaction in metals 

In the framework of many-particle theory of non-transition metals the energy of 
electron subsystem in the field of fixed ions in the adiabatic approximation may be 
written [[|: 

E e = ipo + YsflCR-n) + Y\ X^2(R-n,Rm) + ^ E lf 3 (K n ,R m ,Rl) H (1) 

where each term of the (1) series describes interaction of ion groups through the 
surrounding electron gas and can be represented as a power series in the EIP ||. 
For example 

oo 

^(Ri,R 2 ) = ^^(Rx,^), 



i=2 
oo 



^ 3 (Ri,R 2 ,R3) = E $ 3 ) (Ri,R2,R 3 ), (2) 



i=3 



etc. The indirect interaction of two ions in second order perturbation theory in EIP 
is well known |JTTJ| and equal to 

*?\R) = ^/^ 2r(2 ^)l^)| 2 ^^, (3) 
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where V(q) = — ^jr is the formfactor of the EIP; r^ 2 ^(g) is the sum of two-pole 
diagrams. 

In the third-order perturbation theory in EIP the indirect pair-ion interaction is 
defined by the expression 



3 x 



^\ R ) = 4^4 / d <Hll J J dxV( qi )V(q 2 )V(q 3 )T^( qi , q 2 , q 3 ) 



x 



0-1 



x « (4) 

where T^(qi, g 2 , q 3 ) is the sum of three-pole diagrams; q 3 = (qf + q\ + 2gig 2 a;) 1 / 2 . 

The potential of indirect three-ion interaction after double integration may be 
written [[/J 

oo oo 

<5>f\Ri2, #23, #13) = ^4 / d 9i?i 2 / dq 2 q 2 2 X 



1 1 

x / ^v r (g 1 )y(g 2 )K(g 3 )r (3) (gi,g2,g 3 ) / rf« 



x cos( g( gi^ ^ 12 ^f 2 ! ^ 13 + 92-823*) 
ZK12K23 / , 



x J U1-R12 1 - x 



x Jo (g 2 #23 (l - a: 2 ) V2 (l - z 1 ^ 1/2 ^j , (5) 

where Jq(x) is Bessel function of zeroth order; z = cos(qi,q 2 ); R12, R23, R13- the 
distances between the vertices of a triangle formed by the protons. 

The calculations of potentials were executed at the Wigner-Zeitz radius = 
1,65 which corresponds to zero pressure in the zeroth model of a metal ||. A 
permittivity function in the Heldart-Vosko form was employed. The collections of 
potential values were calculated on set with step equal to 1 Bohr radius for each 
dimension. In the drawings of graphs was used cubic-spline interpolation. 

3. Results 

Pair interionic potential in the third order has a form 

^(R) = e ^ + ^\R) + ^(R) } (6) 

where e is the electron charge. 

Figure 1 shows computed potential <fi*{R) and its components. The curve 1 
corresponds to contribution ^ + $ 2 (R), 2 - $ 2 (R), 3 - <p*(R) ■ Pair interionic 



R 

,(2), 



potential ^ + $ 2 (R) hasn't deep potential well. Stevenson and Ashcroft [|T^] have 
obtained pair interproton potential %■ + $ 2 f° r MH at r$ = 1,6 analogous to 
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showed in Fig. 1 (curve 1). Indirect interionic interaction $2 (i?) has attractive 
nature and forms potential well for nearest neighbouring ion and minimum in the 
repulsive part of the potential <p*(R). The decrease of the density leads to strong 
increase of the deep of minimum in the repulsive part of the potential ip*(R). The 
curve 1 at the Fig. 2 corresponds to potential (p*(R) at r s = 1,72. At r s < 1,65 
the minimum in the repulsive part of tp*{R) turns shallow and its position shifts 
towards smaller R and at rs = 1,55 this minimum disappears (curve 2 in Fig. 2). 
The position of minimum in the repulsive part of <p*{R) in the Fig. 1 corresponds to 
separation R = 1,6ob. It should be mentioned that the separation between nuclei 
in H 2 molecule is equal to 1, 4ag and binding length in H 2 ion is equal to 2a B [T3| . 
The formation of the minimum in the dependence at a distance which is smaller 
than the interproton average distance and strong dependence of the deep of this 
minimum from density of MH may be treated as the tendency to diatomic ordering 
in electron-proton plasma. 

(3) 

The results of calculation of the irreducible three-ion interaction potential $3 
for rs = 1,65 are presented in Fig. 3 in form of the potential relief. The minimal 
interproton separation in metallic phase at P=0 equals to 2, 04a# M. Two ions are 
placed in ordinate axis on distance 2a B and third ion is sited in plane XOY. For 
small distances the potential $ 3 (2, R23, R13) has attractive character. The potential 
$3 3 ' ) (2, R 23 , R13) has anisotropic form and its most rapidly changing takes place in 
OX direction. For large distances the potential $ 3 has oscillatory character. Note 
that $ 3 3) (0,0,0) ~-lRy. 

Fig. 4 (geometry coincides with Fig. 3) shows potential relief for the 3-rd ion 
in the fields of two other ions with pair interactions taken into account, that is the 
potential tp* (R23) + tp* (R13) ■ The minima B and D and a valley which connects them 
correspond to the possible positions of the nearest neighbouring ion in the field of 
two fixed ions (see Fig. 4). It should be noted that ion transition from minimum 
B to minimum D along the valley entails the crossing the potential barrier with 
height ~ 1200 K. This situation remains unchanged when irreducible three-proton 
interaction is taken into account because of the short-acting character of potential 
$ 3 . The quantum consideration of ion motion in the field of two fixed ions results 
to conclusion about splitting of energy level for ion which is the nearest neighbour 
for pair considered. The local minima A and C in Fig. 4 correspond to distances to 
ion 1 and 2 which are less than interionic average distance. 

Fig. 5 shows the potential relief which is created ions pair and taking into account 
pair- and three-ion interactions. It is described by function ip*(R 23 ) + f*{Ri3) + 
$3 3 ' ) (2, R 2 3, -R13). The change of local minima A,B,C,D positions as compared with 

Fig. 4 is unimportant. As a consequence of short-range and anisotropic character of 

(3) 

potential $ 3 the depth of minimum A changes markedly and the depth of minimum 
B increases a bit. The ion wave function in potential relief which is shown in Fig. 
5 has nonzero value near local minima A and C. This means that ions have finite 
probability draw together to separation of the order of internuclei distances in H 2 or 
H 2 ■ From comparison of Fig. 1 and Fig. 5 the conclusion follows that the presence 
of third ion essentially increases this probability. 

Fig. 6 shows computed distribution of conditional probability density for ion 
positions in ions pair field with pair and three-ion interactions taken into account 
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in Boltzmann approximation at temperature T = 1000 K. In this approximation 
distribution of the conditional probability density may be described in the form 



F 1 (R 1 \R 2 ,R 3 )^ex V ^%R 12 ) + i P %R 13 ) + ^\R 12 ,2,R 13 ))/kTj, (7) 

where k is Boltzmann constant. Local minima of the potential relief at Fig. 5 
correspond to peaks of the function Fi(Ri|R 2 , R3). 

It is interesting to consider the potential relief which is created by group of ions. 
Fig. 7 shows potential relief which three ions forming equiangular triangle with side 
2a b and placing in XOY plane create in the same plane towards fourth ion. In 
Fig. 7 the pair-ion interaction is taken into account and potential relief assigns to 
function <£>*(i? 14 ) + ^*{R 2 i) + ^(Rza)- It is easy to see two minima A and B and 
saddle point C in XOY plane. 

Potential relief which is shown in Fig. 8 (the geometry is the same as in Fig. 7 
) and which is described by function 

<P*(Ru) + p*(i2 24 ) + <P*{R 3 i) + $f(2, R 2 4, Ru)+ 

^f ) (R 2 4, 2, #34) + $3 3) (#14, #34, 2) (8) 

takes into account three-ion interactions in the cluster of four ions. Three-ion inter- 
actions transform the saddle C into local minimum and minima A and B become 
deeper considerably. This means that an increase of the particles number in cluster 
lowers the potential barrier to diatomic ordering. 

If three ions are situated so as is shown in Fig. 8 and the fourth ion is placed on 
Z-axis, the dependence of potential with pair- and three-ions interactions taken into 
account is shown in Fig. 9. The presence of deep minimum at z = as is characteristic 
peculiarity of this dependency. Second minimum at z = 3,5cib corresponds to 
equilibrium position of fourth ion. 

The quantum description of ion motion in two-well potential of Fig. 9 yields 
the conclusion about a possibility of tunneling transition from well A to well B. The 
probability of ion transition per 1 s is equal to product of the ion oscillation frequency 
in well A and transmission coefficient through potential barrier. The transmission 
coefficient can be estimated in quasiclassical approach. The ion frequency of zero- 
oscillation in the well A is equal to v = 3, 61 ■ 10 13 Hz. The transmission coefficient 
through potential barrier D = 3, 85 • 10~ 12 . Then ion life time r in well A is equal 
to r = {Dv)~ x = 0, 0072 s. If ion transition into well B takes place the conditions 
for electrons localization on ion pairs will be created. Taking into consideration also 
small characteristic time of an electron subsystem it is possible to consider obtained 
time as characteristic time of the tunneling nucleation of H 2 molecule in MH at 
zero pressure. If the identity ions is taken into account it is possible to make a 
conclusion about manyionic tunneling mechanism for metallic hydrogen transition 
into molecular phase. 



The quantum Monte-Carlo simulations of the molecular hydrogen |L4[ and cal- 
culations of the diatomic phase energy as function of the interproton separation || 
indicate that at chosen density the molecular unit interproton separation is some- 
what less than the neutral molecule separation l,4a#. We calculated the potential 
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relief which is created by ions pair with interproton separation l,5ag (Fig. 10,11). 
Fig. 10 shows potential relief for the third ion in the field of two other ions with 
pair-ion interaction in third order perturbation theory (analogous to Fig. 4). Fig. 
11 (geometry coincides with Fig. 10) shows potential relief for the third ion with 
pair and irreducible three-ion interactions taken into account, that is the potential 
(p*(R 12 ) + p*CR 2 3) + $jj 3) (1.5, #13,^23). As is easily seen from Fig. 10,11 , the ir- 
reducible three-ions interaction essentially deepens the potential well A ( Fig. 11 is 
analogous to Fig. 5 ). The comparison of Fig. 5 and Fig. 11 yields the conclusion 
about strong dependence on the interproton separation i? 12 of the probability of 
triplet ions drawing together to the interproton separation as in molecule H 2 . The 
time of the proton transition from well B into well A (Fig. 11) cannot be considered 
as characteristic time of the tunneling nucleation H 2 molecule in metallic phase since 
interproton separation R\ 2 = l,5ag is less than minimal interproton separation in 
MH phase at P=0 which equals to 2, 04a B 0. 

4. Discussion and conclusions 

Our consideration is based on the assumption that metastable MH phase exists 
at P=0. The calculations show that irreducible three-ion indirect interactions have 
decisive influence on the local structure of the ionic subsystem. The calculations 
of the possible MH structures at P=0 show the tendency for MH to crystallize in 
triangular string with two-dimensional periodicity structure . Our calculations of 

the potential relief which three ions create towards fourth ion confirm this conclusion 

(3) 

and show that taking into account $3 really leads to formation of the hexagonal 

(3) 

structure in the plane of four ions arrangement. Taking into account $ 3 leads to 
transformation of the saddle C (Fig. 7) in local minimum (Fig. 8). Minimum C is 
formed along direction which forms angle 60° with OX axis (the curves on Fig. 7, 8 
are drawn with a step 15° ). 

It is possible to observe an interesting peculiarity of the potential relief which 
ion pair creates towards third ion (Fig. 4,5,6): local minima B and D and valley 
which connects them and has saddle point. To make transition from minimum B to 
D it is necessary to cross over the potential barrier. This means that energy levels 
of the third ion in field of two fixed ions are split. The quantum effect of energy 
levels splitting will be exercised at temperatures which are less than potential barrier 
height. It should be noted here that as a consequence of short- range character of 
potential $ 3 3 ' ) it has weak influence on this effect (Fig. 4,5). It can be supposed 
that the existence of two-level subsystems in amorphous solids at low temperature is 
somehow the display of collective quantum oscillatory states of three particle groups. 

In this study the peculiarities of the pair- and three-ion interactions are inter- 
preted as the vestiges of the molecular phase in metastable MH. These peculiarities 
occur in metastable metallic phase and are absent at megabar pressure (Fig. 2). 
This means that electron-proton system possesses a self-organization effect of the 
"preparing" a transition to the molecular phase. Such effect increases sharply with 
decreasing of the electron gas density. 

The probability of ions drawing together to the interproton separation as in 
molecule H 2 depends essentially on the electron gas density. With a density de- 
creasing the probability of the manyionic tunneling nucleation of the molecules H 2 
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sharply increases. Therefore some interest has a study of the life time of metastable 
MH relative to the homogeneous nucleation of the molecular phase in pressure inter- 
val < P < Pt (Pt is the pressure of the molecular hydrogen transition in metallic 
phase). 

The calculation of the life time of the MH as macroscopic system needs the con- 
sideration both the mechanism of homogeneous nucleation of molecular phase and 
kinetics of formation of the insulating phase droplets in the metastable metal. In 
this paper we did not consider a kinetic stage. Such consideration must include a 
search of volume and surface free energy minimum of heterogeneous system consist- 
ing of metastable metal and insulating droplets. In accordance with the classical 
and quantum theories of the new phase formation . ]I5| , [16j calculated in this 



paper life time of metallic phase relatively tunneling nucleation of the H 2 molecule 
may differ from life time of MH as macroscopic system by a factor ~ 3-5. The new 
experimental data on structures and phase diagrams for molecular hydrogen and 
deuterium at megabar pressure |T7J, |TB[ may be very useful by consideration of the 



kinetic stage of the formation and growth of the insulating phase droplets in the 
metastable MH. 
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Figure captions 

Fig. 1. Potential <p*(R) and its components for r s = 1,65. 

Fig. 2. Potential (p*(R) for r s = 1, 72 and r s = 1, 55. 

Fig. 3. Irreducible three-ion interaction potential $3 (2, i?i3, i?23)- 

Fig. 4. Potential relief for the 3-rd ion in the field of two ions with pair interac- 
tions taken into account. 

Fig. 5. Potential relief for the 3-rd ion in the field of two ions with pair- and 
three-ion interactions taken into account. 

Fig. 6. Distribution of conditional probability density for ion positions in ions 
pair field. 

Fig. 7. Potential relief which three ions create toward fourth ion (pair interactions 
taken into account). 

Fig. 8. Potential relief which three ions create toward fourth ion (pair- and 
three-ion interactions taken into account). 
Fig. 9. The dependence of potential $(z). 

Fig. 10. Potential relief for the 3-rd ion in the field of two ions (interionic 
separation is equal to l,5a B ) with pair interactions taken into account. 

Fig. 11. Potential relief for the 3-rd ion in the field of two ions (interionic sepa- 
ration is equal to 1,5ob) with pair- and three-ion interactions taken into account. 
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